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ABSTRACT 


Spatial  correlations  of  amplitude  anomalies  have 
been  conducted  over  LASA  and  LASA  subarrays  to  test 
the  hypothesis  that  these  anomalies  exhibit  spatial 
stationarity .  The  evidence  indicates  that  the  anomaly 
process  cannot  be  considered  to  be  covariance  stationary 


I  o  INTRODUCTION 


This  report  describes  results  of  part  of  a  study  of  ampli¬ 
tude  anomalies  at  LASA.  The  object  of  this  part  of  the  study  was 
to  see  if  these  anomalies  could  be  treated,  in  some  sense,  as  a 
stochastic  process . 

(1  2 ) 

It  has  been  demonstrated  '  that  the  normalized  short  period 
peak-to-peak  amplitudes  of  teleseismic  events  have  a  log  normal 
distribution.  That  is,  if  the  amplitudes  of  a  set  of  events  from 
the  same  geographic  region  are  measured  and  their  logarithms  are 
taken,  ther  we  find  that 


log  a . .  =  log  L . .  -  - 
13  13  N 


N 

E 

-l-l 


log  L 

i] 


In  this  equation,  the  L, .  are  either 

1D 


has  a  normal  distribution, 
the  measured  peak-to-peak  amplitudes  at  all  elements  in  a  LASA  sub¬ 
array  or  are  the  peak-to-peak  amplitudes  observed  at  the  center 
elements  of  the  subarrays.  The  index  j  is  on  the  seismometer  and 
i  is  an  event  index. 

The  distribution  of  log  amplitudes  is  not  normal  if  the  col¬ 


lection  includes  all  elements  in  LASA. 


The  variance  of  log  a. .  is 


larger  in  the  case  where  the  L^_,  are  the  observed  amplitudes  at 
the  center  elements.  The  variances  are  the  same  at  each  subarray 
when  the  are  the  amplitudes  of  the  elements  in  a  subarray.  Thus, 
these  variances  can  be  pooled  after  normalization. 

These  anomalies  are  assumed  to  be  real  in  that  a  precisely 
repeated  event  should  produce  the  same  amplitudes  at  the  seis¬ 
mometers  as  the  original.  The  anomalies  vary  however  for  events 
from  the  same  geographic  region  and  it  is  unlikely  that  a  calibra¬ 
tion  of  the  earth  would  be  a  practical  procedure  with  which  to 
eliminate  anomaly  effects.  Rather,  a  statistical  approach  may  be 
a  more  reasonable  way  to  proceed. 
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The  fact  that  the  anomalies  in  the  subarrays  can  be  pooled 
after  normalization  suggests  that  one  may  successfully  hypothesize 
that  these  anomalies  exhibit  Sx^atial  stationarity<>  That  is,  although 
there  may  be  slowly  varying  anomaly  effects  with  distance*,  with 
these  removed  the  expectation  of  a  particular  amplitude  anomaly  is 
independent  of  spatial  location. 

Beyond  this  we  desire  that  the  anomaly  process  be  covariance 
stationary.  If  this  is  so  th'in  the  covariance  function  will  serve 
as  a  measure  of  the  distance  which  should  be  placed  between  seis¬ 
mometers  so  that  they  will  exhibit  independent  amplitude  estimates. 
Further,  since  the  anomalies  are  log-normal  no  other  statistic  is 
needed  since  the  covariance  function  is  a  complete  statistic  for 
normally  distributed  variables . 

2.  SPATIAL  CORRELATIONS  OVER  LASA 

This  section  investigates  the  possibility  of  correlation 

among  the  peak-to-peak  amplitudes  across  all  of  LASA.  Data  used 

were  from  eleven  Fiji  Island  earthquakes  which  occurred  at  243° 

azimuth  and  from  9,500  km  to  10,500  km  distance  from  the  center  of 
(3) 

LASA.  From  these  eleven  events,  correlation  coefficients  were 

computed  as  spatial  displacements  were  made  over  LASA.  In  the 
computations,  the  logarithms  of  the  normalized  data  are  used  for 
the  reasons  cited  in  the  introduction  of  this  report.  We  define 
the  estimate  of  the  coefficient  of  correlation  to  be 

r  =  (£  X.Y.)  /  (I  Xi2)Js  (E  (1) 


*  We  note,  however,  that  on  a  larger  scale,  these  slow  variations 
also  are  distributed  as  log  normal. 
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r  ...  *T 


i.  e.  we  have  set  the  mean  of  X  and  of  Y  to  be  zero.  The  justifi¬ 
cation  for  pre-determining  X  and  Y  is  as  follows  %  the  estimate  of 
the  coefficient  of  correlation  is  generally  defined  to  be 


A 

Cov  (X,  Y) 


which  for  small  samples  turns  out  to  be 


(N-i)  E  (X..-X)  (Yi-Y) 


E(Xi-X)  (Yi-Y) 

1  r- - 

Vz(X  -X)2  a/e^-Y)2 


Now  r  is  the  coefficient  of  correlation  for  the  two  dimensional 
distribution  of  (X,  Y) .  But  the  instruments  that  make  up  X  and  Y 
are  actually  subsets  of  the  set,  say,  Z  which  contains  the  21  Bub- 
array  center  seismometers .  Since  X  and  Y  are  samples  of  Z,  and 
since  we  seek  the  mean  of  the  population  from  which  X  and  Y  were 
drawn,  it  is  reasonable  to  use  Z  as  a  better  estimate  of  the  pop¬ 
ulation  mean  than  that  given  by  either  X  or  Y.  Due  to  the  normal¬ 
izing  process,  Z  is  zero.  Setting  X  =  Y  =  Z  =  0  in  (2)  yields  the 
equation  given  by  (1).  In  addition,  another  degree  of  freedom  is 
gained  in  estimating  the  coefficient  of  correlation  since  the  mean 
is  not  estimated.  We  will  apply  equation  (1)  throughout  the  re¬ 
mainder  of  this  report  for  all  calculations  of  the  coefficient  of 
correlation. 


A.  LASA  Spatial  Correlations  Using  Average  Log  of  Normalized 
Peak-To-Peak  Amplitudes . 

Since  all  eleven  events  were  relatively  closely  grouped  in 
comparison  to  the  overall  path,  the  average  of  the  logarithms  tof 
the  normalized  peak-t-  ^-peak  amplitudes  at  a  given  subarray  was 
used  as  an  estimate  of  the  true  value  for  that  subarray  (cf  .Table  I) 


erJ&i, 


—SB*. 


These  average  values  were  used  in  the  calculations  of  the 
spatial  correlation  coef f icients „  Correlations  were  computed 
along  a  line  parallel  to  the  incoming  signal  (243°  az.)  and 
another  set  along  a  line  perpendicular  to  the  incoming  signal 
(153°  az.).  Due  to  the  configuration  of  LASA,  few  displacements 
exist  where  enough  subarrays  intersect  to  compute  a  valid  coef¬ 
ficient  of  correlation.  Table  II  lists  the  subarrays  which  were 
selected  as  intersecting  at  various  displacements.  Figure  1 
presents  the  coefficients  of  correlation  plotted  against  the 
spatial  shifts  using  the  average  logarithm  of  the  normalized 
peak-to-peak  amplitudes  over  all  events  as  the  estimate  of  the 
true  value  for  each  subarray.  The  numerical  values  are  shown  in 
Table  III,  Part  A. 

B .  IAS A  Spatial  Correlations  of  Single  Events 

The  correlation  coefficients  for  the  individual  events  were 
computed.  That  is,  each  event  was  consx  'red  singly  as  the  spatial 
shifts  were  made  across  LASA  at  153°  and  243°  azimuth  as  in  A. 

The  coefficients  are  presented  in  Table  III,  Part  B.  Figures  2 
through  11  show  the  behavior  of  the  coefficients  as  the  spatial 
displacements  are  made.  Note  that  Table  III,  Part  C,  contains 
the  average  coefficient  of  correlation,  r,  for  various  displace¬ 
ments,  and  Figure  12  presents  the  graph  of  r  vs.  the  spatial 
shifts.  It  is  of  interest  to  compare  the  graphs  of  Figures  1 
and  12,  i.e.  the  coefficient  of  correlation  of  the  average  normal¬ 
ized  amplitudes  vs.  displacement  against  the  average  coefficient 
of  correlation  of  the  individual  events  vs.  displacement. 

3,  SPATIAL  CORRELATIONS  OVER  LASA  SUBARRAYS 

After  examining  correlations  over  all  LASA,  we  directed  our 
attention  to  spatial  correlations  over  subarrays.  The  concept  of 
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spatial  shifting  over  a  subarray  is  the  same  as  that  over  all  of 
LASA  except  that  unlike  LASA  the  configuration  of  a  subarray  is 
well  defined  in  terms  of  concentric  circles  and  radial  spokes 
every  60  .  The  spatial  shifts  over  the  subarrays  were  made  along 
the  legs  (radial  spokes)  in  this  instance  rather  than  with  respect 
to  the  origin  of  the  event . 

A.  Subarray  Spatial  Correlations  for  Single  Events 

Three  N.  Colombia  events  (21  Dec.  65?  21  April  66?  12  June 

66)  were  chosen  on  the  basis  of  availability  of  complete  tapes  of 
all  525  instruments  in  LASA  for  events  which  originated  very  cisse 
together.  Again  the  logarithms  of  the  normalized  data  were  used 
to  compute  the  coefficient  of  correlation  and  the  means  of  the 
data  samples  were  set  to  zero  as  motivated  in  Section  2,  LASA. 

As  spatial  displacements  were  made  across  each  subarray  (at  0.5  km 
increments),  coefficients  of  correlation  were  computed  and  graphs 
of  coefficients  vs.  displacement  were  drawn  for  each  leg  for  every 
subarray.  Some  representative  graphs  are  given  in  Figures  13 
through  18.  The  patterns  in  the  graphs  were  cross-checked  among 
the  same  subarrays  for  all  three  events  and  among  similarly  oriented 
subarrays  for  the  same  event.  No  consistent  relationships  were 
discovered  with  this  approach.  Contouring  (Figures  19  and  20  are 
two  examples)  seemed  to  suggest  certain  patterns  (viz.,  that  the 
contour  "pointed"  towards  the  direction  of  the  event)  and  some 
consistencies  were  found,  but  the  presence  of  exceptional  and 
contradictory  contours  made  such  a  conclusion  at  best  doubtful. 

It  was  thought  that  teleseismic  explosions  occurring  close  together 
might  be  better  sources  to  study  for  contouring  in  this  manner. 

This  step  is  taken  in  3-C  below,  but  first  we  will  mention  some 
additional  investigation  of  the  three  N.  Colombia  events. 
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B-  Testing  For  Correlation  Among  Three  Colombia  Earthquakes 
Coefficients  of  correlation  were  computed  over  each  subarray 
for  paired  events  using  the  normalization  described  earlier  and 
equation  (1)„  That  is,  the  twenty-five  readings  of  the  A  subarray 
for,  say,  the  21  December  1965  earthquake,  were  tested  against  the 
respective  twenty-five  readings  cf  this  subarray  for,  say  the  21 
^Pril  1966  earthquake „  The  three  independent  earthquakes  yield 
three  distinct  pairs  for  testing  in  this  manner,  and  each  pair 
of  events  has  a  maximum  of  twenty-one  coefficients  to  be  computed... 
one  coefficient  corresponding  to  each  subarray. 

In  nearly  every  case,  the  coefficient  was  significant  (i.e. 
correlation  existed)  and  in  many  cases  (nearly  two-thirds)  the 
coefficients  were  between  .8  and  1.0.  Table  IV  summarizes  the 
results  of  these  calculations.  This  test  can  be  considered  as 
examining  the  correlation  among  the  instrument  responses  .3^  the 
source  distance  is  varied.  In  the  next  step,  the  converse  of 
this  procedure  was  done  . . .  the  respective  elements  for  similarly 
oriented  subarrays  (account  being  taken  for  long  and  short  con¬ 
figurations)  were  inspected  for  correlation  for  each  earthquake. 

This  time,  correlations  were  not  detected  with  the  exception  of 
one  case  which  may  be  ascribed  to  chance  (cf.  Table  V). 

c  0  Subarray  Spatial  Correlations  For  Two  Shots 
Two  nuclear  shots  from  Kazakh  (November  21,  1965;  February  13, 
1966)  were  selected  to  compare  with  the  results  of  the  three 
Colombian  earthquakes.  The  correlation  coefficients  were  contoured 
as  had  been  done  of  the  earthquakes,  and  a  search  was  made  to  detect 
similar  patterns  between  the  same  subarrays  for  the  two  shots.  The 
number  of  similarities  found  was  no  greater  than  that  expected  purely 
by  chance. 


D .  Testing  For  Correlation  Between  Two  Kazakh  Shots 


The  coefficient  was  computed  over  each  subarray  for  both 
events  as  was  described  in  3-B.  Six  of  the  subarrays  registered 
a  positive  correlation,  one  a  negative,  and  the  remaining  thir¬ 
teen  (one  subarray  was  inoperative)  no  correlation,  as  is  shown 
in  Table  VI.  Performing  the  converse  test  on  this  data  provided 
results  similar  to  those  obtained  from  the  earthquakes,  i.e.,  no 
correlations  were  detected  among  similarly  oriented  subarrays 
for  the  same  event  (cf  Table  VII ) . 

4.  CONCLUSIONS 

Only  tentative  conclusions  can  be  drawn  from  this  data.  The 
sparce  sampling  of  the  LASA  array  limits  the  reliability  of  the 
correlation  coefficients  which  were  computed.  For  this  reason  a 
uniformly  spaced  grid  of  seismometers  would  have  aided  this  study. 

It  is  likely  that  the  anomaly  process  cannot  be  considered  to 
be  spatial  covariance  stationary.  Since  this  process,  is  in  fact, 
a  description  of  the  underlying  gsology  one  might  have  hypothe¬ 
sized  this  from  the  beginning.  The  author  does  not  have  a  simple 
explanation  for  the  log-normal  distribution  of  the  anomalies  or 
for  their  apparent  stationarity  although  this  too  reflects  the 
geology. 
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TABLE  II 


Intersecting  Subarrays  For  Various  Displacement  Dis¬ 
tances  and  Azimuths  Used  In  Spatial  Correlation 


Direction  of  Displacement 
+  10  km 


Cl,  B4 
Bl,  AO 
C2,  B2 
B2 ,  C3 
AO,  B3 
B4 ,  C4 


Direction  of  Displacement 


+  10  km 

B3, 

C4 

C3, 

B3 

C2 , 

Bl 

AO, 

B4 

Bl, 

Cl 

+  80  km 

E3, 

E4 

F2 , 

D2 

B4 , 

F4 

-  24'° 


+  18  km 

+  22  km 

B3,  D3 

AO,  D3 

Bl,  B3 

Cl,  C4 

Dl,  Bl 

Bl,  B3 

Cl,  C4 

C2,  C3 

C2 ,  C 3 

Dl,  Bl 

153° 

+  17  km 

+  km 

D2 ,  B2 

D2,  B4 

B2,  Bl 

B2,  D4 

C  3 ,  C4 

B4 ,  D4 

C2  ,  Cl 

Dl,  El 

+  90  km 

F4 ,  AO 
El,  F2 
C3,  E2 


This  table  lists  the  intersecting  subarrays 
for  various  displacement  distances  and  azimuths 
which  were  used  in  the  spatial  correlation 
described  in  the  report  text. 


TABLE  III 


Coefficients  of  Correlation  For  Individual 
Events  At  Designated  Spatial  Shifts 


PART  B 


Event  Nos. 

152 

-.403 

-.004 

-.112 

.341 

-.401 

X 

-.950 

X 

171 

-.121 

-.207 

-.221 

.108 

-.829 

X 

X 

-.826 

197 

-.877 

.057 

.141 

.469 

.146 

X 

X 

X 

219 

-.389 

-.724 

-.753 

.374 

.381 

X 

X 

X 

238 

-.029 

.224 

.395 

-.035 

-.530 

-.979 

X 

X 

247 

.491 

.236 

.236 

.194 

-.689 

-.567 

X 

-.330 

253 

.592 

.894 

.894 

.227 

-.797 

-.273 

-.773 

X 

254 

-.206 

-.615 

-.615 

.328 

-.506 

X 

X 

X 

271 

.880 

.820 

.820 

.085 

-.921 

x 

X 

.727 

359 

-.223 

-.686 

-.686 

-.066 

-.828  j 
_ 1 

X 

X 

-.946 

PART  C 


Avg.  Coef. 
of  Correl. 

-.028 

-.001 

-.010 

1 

.203 

-.498 

-.606 

-.861 

-.344 

cf  .Fig. 

of  Each 
Event,  r 

1 

12 

. . 


j  TABLE  IV 

Subarray  correlations  for  the  Colombia  Earthquakes 


21  Dec  65 

21  Apr  66 

21  Dec  65 

12  Jun  66 

21  Apr  66 

12  Jun  66 

Computed 
__E _ 

Critical 

Computed 

r 

Critical 

r,5* 

Computed 

r 

Critical 
r ,  5% 

B1 

0.693 

0.388 

0.613 

0.388 

0.537 

0.388 

r3 

0.955 

0.388 

0.851 

0.388 

0.776 

0.388 

T4 

0.987 

0.388 

-0.477 

0.388 

-0.503 

0.388 

AO 

o.4ie 

0.388 

0.637 

0.388 

0.235 

0.388 

83 

0.815 

0.388 

0.754 

0.308 

0.959 

0.388 

C4 

0.733 

0.396 

0.933 

0.388 

0.943 

0.396 

84 

0.81S 

0.388 

0.754 

0.396 

0.790 

0.396 

Cl 

0.939 

0.‘388 

0.893 

0.388 

•  1 

0.844 

0.388 

C2 

0.898 

0.388 

0.782 

0.388 

0.639 

0.388 

82 

0.857 

0.388 

0.767 

0.388 

0.945 

0.388 

C3 

0.902 

0.388 

0.746 

0.388 

0.836 

0.388 

03 

0.986 

0.388 

0.977 

0.388 

0.979 

0.388 

04 

0.563 

0.388 

0.368 

0.388 

0.634 

0.388 

Dl 

0.831 

0.404 

0.882 

0.404 

0.931 

0.388 

02 

0.412 

0.396 

0.452 

0.388 

0.511 

0.396 

E3 

0.945 

0.388 

0.815 

0 . 388 

0.843 

0.388 

E4 

0.450 

0.388 

0.511 

0.388 

0.805 

0.388 

El 

0.900 

0.388 

0.825 

0.388 

0.806 

0.388 

FI 

0.956 

0.388 

0.916 

0.388 

0.896 

0.388 

E2 

0.873 

0.388 

0.889 

0.388 

0.901 

0.388 

F2 

0.956 

0.388 

0.860 

0.388 

0.860 

0.388 

Subarray  correlations  for  the  Colombia  earthquakes 


The  critical  values  are  determined  by  using  the 
distribution  adiere 


A  detailed  explanation  is  presented  in  Snedecor's 
"Statistical  Methods"  Fifth  Edition,  pp  173,  174. 


TABLE  V 


TABLE  VI 


Subarray  Correlations  For  The  Two  Kazakh 


Subarray 

Designation 

B1 

F3 

F4 

AO 

B3 

C4 

Cl 

C2 

B2 

C  3 

D3 

D4 

D1 

D2 

E3 

E4 

El 

FI 

E2 

F2 


Coeff.  of 
Correlation 

0.067 

0.021 

-0.042 

0.615 

0.898 

0.905 

0.276 

-0.264 

-0.157 

-0.030 

0.275 

0.300 

0.333 

0.071 

0.554 

0.788 

0.044 

-0.573 

0.595 

0.220 


Explosions 

Critical 
r ,  5% 

0.388 

0.388 

0.388 

0.388 

0.388 

0.388 

0.388 

0.388 

0.388 

0.3F-6 

0.388 

0.388 

0.388 

0.388 

0.388 

0.388 

0.388 

0.398 

0.388 

0.388 


Subarray  correlations  for  the  two  Kazakh  explosions 


TABLE  VII 

Correlations  of  Similarly  Oriented  Subarrays 
For  The  Kazakh  Explosions 


21  Nov. 

1965 

13  Feb. 

1966 

Computed 

r 

Critical 
r,  5% 

Computed 

r 

Critical 
r ,  5% 

C4,  El 

-0.075 

0.388 

-0.062 

0.388 

Bl,  C2 

0.068 

0.388 

0.021 

0.388 

Bl.,  D2 

0.159 

0.388 

-0.282 

0.  388 

Bl,  E3 

0.170 

0.388 

0.024 

0.388 

C2,  D2 

0.229 

0.388 

0.294 

0.388 

C2 ,  E3 

0.214 

0.388 

-0.175 

0.383 

02,  E3 

-0.240 

0. 388 

-0.117 

0.388 

nti  PMALici  to  oinernoii  or  cvorr 


O  displacements  parallel  to  direction  or  EVENT . 


Figure  2.  Event  152 


DISPLACEMENTS  PARALLEL  TO  DIRECTION  OP  EVENT . 


.  f 


O  displacements  parallel  to  direction  of  event  . 


Figure  4.  Event  197 


DISPLACEMENTS  PARALLEL  TO  DIRECTION  OF  EVENT  , 


DISPLACEMENTS  parallel  to  direction  or  EVENT. 


displacements  parallel  to  direction  of  event  . 


Figure  7.  Event  247 


O  displacements  parallel  to  direction  of  event. 


Figure  8.  Event  253 


s _ *  •  '  *  ■  •  ^  •'  •  \  .  <•*$ 


Figure  11.  Event  359 
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Figure  13.  N.  Colombia,  21  December  x965 
Seismogram  No.  6393, C4, Leg  1 


Figure  14.  N.  Colombia,  21  December  1965 
Seismogram  No.  6393, C4, Leg  2 
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N.  Colombia,  12  June  1966 
Seismogram  No.  7635,  C4, 
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Figure  18.  N.  Colombia,  12  June  1966 
Seismogram  No.  7635,  C4, 
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